High-throughput electronic band structure calculations: challenges and tools 



Wahyu Setyawan 1 and Stefano Curtarolo 1,2 

Dept. of Mechanical Engineering and Materials Science and Dept. of Physics, Duke University, Durham, NC 27708. 

2 corresponding author: stefano@duke.edu 



Abstract 

The article is devoted to the discussion of the high-throughput approach to band structures calculations. We present 
scientific and computational challenges as well as solutions relying on the developed framework (Automatic Flow, 
AFLOW/ACONVASP). The key factors of the method are the standardization and the robustness of the procedures. 
Two scenarios are relevant: 1) independent users generating databases in their own computational systems (off-line 
approach) and 2) teamed users sharing computational information based on a common ground (on-line approach). 
Both cases are integrated in the framework: for off-line approaches, the standardization is automatic and fully inte- 
grated for the 14 Bravais lattices, the primitive and conventional unit cells, and the coordinates of the high symmetry 
k-path in the Brillouin zones. For on-line tasks, the framework offers an expandable web interface where the user 
can prepare and set up calculations following the proposed standard. Few examples of band structures are included. 
LSDA+U parameters (U, J) are also presented for Nd, Sm, and Eu. 

Keywords: High-Throughput, Combinatorial Materials Science, Computer Simulation, Brillouin Zone Integration, 
VASP, AFLOW, ACONVASP 



Over the past decade computational materials sci- 
ence has undergone a tremendous growth thanks to 
the availability, power and relatively limited cost of 
high-performance computational equipment. The high- 
throughput (HT) method, started from the seminal pa- 
per by Xiang et ah for combinatorial discovery of su- 
perconductors 1 1], has become an effective and efficient 
tool for materials development J2, S B B 0] and pre- 
diction A S, B [H El E Q. Recent examples of 
computational HT are the Pareto-optimal search for al- 
loys and catalysts HQ, the data -mining of quantum 
calculations method leading to the principle-component 
analysis of the formation energies of many alloys in 
several configurations H3, ITTl Tl^L [l6[ \vh , the high- 
throughput Kohn-anomalies search in ternary lithium- 
borides 1181 [T9L I20IL and the multi-optimization tech- 
niques used for the study of high-temperature reactions 
in multicomponent hydrides 112 l[|22l 12311. 

In its practical implementation, HT uses some sort 
of automatic optimization technique to screen through 
a library of candidate compounds and to direct further 
refinements. The library can be a set of alloy prototypes 
ll24l[l2ll or a database of compounds such as the Pauling 
File IS or the ICSD Database 00. An important 
difference between the several-calculations and the HT 



philosophies is that the former concentrates on the cal- 
culation of a particular property, while the latter focuses 
on the extraction of property correlations which are used 
to guide the search for systems with ad-hoc character- 
istics. The power of HT comes with a cost. Due to the 
enormous amount of information produced, standard- 
ization and robustness of the procedures are necessary. 
This is especially true if one were concomitantly opti- 
mizing thermodynamics and electronic structure, which 
is required, for instance, in catalyst design 12811 . in ac- 
celerated "battery materials" discovery Q29JL and super- 
conducting materials development tl9l|20|]. Therefore 
a rational HT computational framework must contain a 
general, reliable, and standardized electronic structure 
analysis feature. It must determine the symmetry auto- 
matically, the Brillouin zone (BZ) integration path for 
all the possible 14 Bravais Lattices with all their vari- 
ous sub-cases, and put the direct and reciprocal lattice 
vectors in the appropriate standardized form, so that 
data can be exchanged and recycled between different 
projects. Although Brillouin zones integration paths 
have been included in books and literature for the last 
few decades [3^ Ell Ell 0L a standardized defini- 



tion of the paths for all the different cases is, to the best 
of our knowledge, missing. 
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In this article, we describe the BZ paths features 
of AFLOW 113 511 . which is our free framework for 
performing high-throughput thermodynamics and elec- 
tronic structure calculations on top of DFT ab initio 
codes (currently the Vienna Ab-initio Simulation Pack- 
age (VASP) but the porting to other DFT packages, such 
as Quantum Espresso 113611 is underway). A typical task 
involves structural optimizations to a ground state ("re- 
laxation run"), determination of charge density and its 
projection onto electronic orbitals ("static run"), and 
calculation of energy levels along a path of "important" 
wave vectors ("bands run"). We refer this set of k-points 
as k-path. 

Some definitions are pertinent for automatic con- 
struction of the k-path. A k-point is a symmetry point 
if its site symmetry contains at least one point symme- 
try operation that does not belong to the site symmetry 
of the neighboring points in sufficiently small vicinity. 
Similarly, a line or a plane forms a symmetry line or 
a symmetry plane if it contains at least one symmetry 
point and all of the k-points on the line or plane have site 
symmetry with at least one point symmetry operation 
not possessed by the neighboring points. The k-path 
must be carefully chosen so that the electronic prop- 
erties of a materials imposed by its underlying crystal 
symmetry are correctly described. For example, GeF4 
(cIlO, ICSD #202558, space group #217, 143m) has an 
indirect gap. The conduction band minimum (CBM) 
occurs at point T, while the valence band maximum 
(VBM) occurs at the point H of the BZ for the body- 
centered cubic (BCC) lattice (Figure [28]). If a cubic unit 
cell were used instead the BCC, one would incorrectly 
find a direct gap occurring at point Y. 

The coordinates of symmetry k-points are more con- 
veniently expressed as fractions of the reciprocal lattice 
vectors. Therefore, the primitive lattice vectors needs 
to be properly defined in a standardized fashion. In or- 
der to build a primitive cell, the framework AFLOW per- 
forms the following procedure. 

i. Given any input structure (unit cell and the coordi- 
nates of the basis atoms), AFLOW reduces it into a min- 
imal set of basis atoms in a primitive cell. 

ii. A set of symmetry properties are calculated: lattice 
point group, crystal point group, crystal family, factor 
group, space group operations, Pearson symbol, and the 
Bravais lattice type. If the input structure contains in- 
formation about the space group number, it will be used 
to double check the Bravais lattice type. 

iii. A standard conventional cell is then identified and 
constructed, and whenever possible, the lattice vectors 
are ordered according to the axial lengths and the inter- 
axial angles. This ordering eliminates some choices in 



the possible shapes of BZ in certain Bravais lattices. For 
example, let us consider the body-centered orthorhom- 
bic lattice. Depending on the ratios of the axial lengths 
the body-centered orthorhombic cell has three possible 
shapes of the BZ, and consequently three different co- 
ordinates of the symmetry points. By ordering the con- 
ventional lattice vectors so that |ai| < \^\ < \^\ (i.e. 
a < b < c) a unique choice remains. Furthermore, the 
ordering improves the similarity between band struc- 
tures of different compounds with the same lattice: the 
proportion between the length of each path segment in 
the band structure will be comparable, 
iv. From the ordered conventional unit cell, a standard 
primitive cell is created. Amongst all the possible prim- 
itive cells, AFLOW, chooses the one with the reciprocal 
lattice vectors passing through the center of the Bragg 
planes belonging to the first BZ. The choice has consid- 
erable practical advantages. It enforces the reciprocal 
lattice vectors to be as perpendicular as possible within 
each other (Minkowski lattice reduction | 37]) and min- 
imizes the number of the plane waves basis set used in 
the quantum mechanical code, allowing faster conver- 
gence and smaller memory requirements. Q After all 
these steps, the standard primitive unit cells of the 14 
Bravais lattices as calculated by the AFLOW package are 
safe to be used in the "relaxation" and "static" calcula- 
tions. As usual, the latter calculation will be used for 
the electronic density of states. 

The k-paths used in the band structure analysis are 
constructed from the irreducible part of the first Bril- 
louin zone (IRBZ). A symmetry line is included in the 
path if it belongs to the edges of the IRBZ, otherwise 
it is included only if it carries one or more new point 
symmetry operations with respect to those of its ex- 
tremes. Duplicate lines, due to the reciprocal lattice 
point group symmetry and translations are also omit- 
ted. A point possessing only identity (E) and inversion 
(/) operations can not form a symmetry line, however, it 
may still be a zero-slope point H38M . Examples are point 
L in face-centered orthorhombic, points N, Ni , and M 
in C-centered monoclinic, and all symmetry points in 
triclinic. To illustrate, TIF (oF8, ICSD #30268, space 
group #69, Fmmm) has CBM and VBM occuring at 
point L (Figure [34]). This result would have not been 
obtained if point L were excluded. Therefore, for com- 
pleteness, a line from point T to such point is included 
in the k-path. 



lr The common visualization softwares used for generating BZs 
simply take the input lattice vectors, generate reciprocal vectors and 
produce BZs which are not necessarily the Wigner-Seitz cells of the 
reciprocal lattice. 
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Notably for triclinic, monoclinic, and rhombohedral 
systems, the shape of Wigner-Seitz cell of the reciprocal 
lattice depends nontrivially on the lattice vectors. For 
this reason, some researchers use the parallelepiped of 
primitive reciprocal lattice vectors, centered at k = 0, to 
define a plausible BZ. Even though the energy is con- 
tinuous throughout such BZ l30ll . its faces, in general, 
are not parallel to the symmetry planes of the lattice. 
Therefore, there is no guarantee that any line connect- 
ing two k-points on a BZ face will be a symmetry line. 
In addition, since parallelepiped unit cells have only 8 
points, one at each corner of the IRBZ, one would miss 
some symmetry points on the Wigner-Seitz cell of re- 
ciprocal lattice. Consequently, a complete irreducible 
set of symmetry lines would not be obtained. For ex- 
ample if we were using a parallelepiped as BZ in the C- 
centered monoclinic variation MCLCi (Figure [T7]), the 
site symmetries of points l h^i, V2D2, l h(^i + D3), and 
V2 (t>2 + l>3) would be only E and /. Therefore, the line 
V2bi-Y would not be, for example, a symmetry line. 
Furthermore, points like X, Xi, I, and Ii and their re- 
lated symmetry lines (C2 about x-axis) would not be 
included in such simplification. To conclude, although 
we believe that the simple parallelepiped BZ can be use- 
ful in some particular difficult cases, we think that the 
solution is not appropriate for a full automatic and high- 
throughput implementation of the band structure analy- 
sis. For this reason, inside the framework AFLOW, all 
the BZ and their k-paths are derived from the Wigner- 
Seitz cell of the reciprocal lattice and the available sym- 
metries. In Appendix, for all the Bravais lattices and 
variations, we present the conventional and primitive 
lattice vectors implemented in AFLOW, the coordinates 
of high- symmetry k-points for the path, the shape of the 
BZ, and an example of band structure calculation for a 
compound extracted from the ICSD database. 
Off-line implementation. 

The effort in developing the standardized tool 
AFLOW comes from the ongoing generation of an ex- 
tensive database of electronic band structure for in- 
organic crystals for scintillator materials design [39]. 
We have extracted approximately 195,000 structures 
from the Inorganic Crystal Structure Database (ICSD) 
E El HE H. AFLOW is equipped with utilities to 
select structures of interest. Selection criteria can be 
based on atomic number or element's name, mass den- 
sity, number of atoms per primitive unit cell, chem- 
ical formula, structure prototype, space group num- 
ber, lattice type, ICSD entry number, etc. Features to 
remove/include structures containing certain elements, 
partial occupancies, and redundancies (structures with 
the same chemical formula and space group number) are 



also implemented. Each structure is given a label which 
is composed of the structure's chemical formula (in al- 
phabetic order) and the ICSD entry number. This label 
is the only information that AFLOW requires to produce 
the band structure. After the structure selection is per- 
formed, a list of labels is produced. Based on such la- 
bels, AFLOW creates a subdirectories for each structure 
and the necessary input file for the band structure calcu- 
lation with VASP (porting to other DFT packages, such 
as Quantum Espresso B36I1 is underway). For running 
the DFT package, AFLOW has an option to run only 
one structure and exit, or to search through subfolders 
and run those that have not been calculated yet, or to 
wait for new structures to run. If started as a common 
Unix daemon through the queue of a computer cluster, 
AFLOW will generate, run, correct and converge many 
calculations per day, with minimum human input. 
Electronic structure database implementation. 

The database of band structures under construction 
is calculated using VASP within the General Gradient 
Approximation of the density functional theory 14211 . 
We use projector augmented waves pseudopotentials 
with exchange correlation functionals as parameterized 
by Perdew-Burke-Ernzerhof |43[ Hi . All structures 
are fully relaxed with a convergence tolerance of 1 
meV/atom using dense grids of 3,000-4,000 k-points 
per reciprocal atom for the integrations over BZ. A 
much denser grid of 10,000 is implemented for the 
static run to get accurate charge densities and density 
of states. Monkhorst-Pack scheme [45] is employed 
in the grid generation except for hexagonal and rhom- 
bohedral systems in which T-centered grid is selected 
for faster convergence. At the beginning of relax- 
ation, a spin-polarized calculation is performed for all 
structures. Then, if the magnetization is smaller than 
0.025 jug/atom, the spin is turned off for the next re- 
laxations and subsequent calculations to enhance the 
calculation speed. At the completion of each calcula- 
tion ("relax"^"relax"^"static"— >"bands"'), appropri- 
ate MATLAB scripts are invoked for data analysis and vi- 
sualization. All these steps are done automatically. One 
of the most difficult challenges in the high-throughput 
combinatorial search is about the response to erroneous 
interruption of the one or more of the flows work- 
ing on a big set of problems, concurrently. The most 
common cause is insufficient hardware resources. Pre- 
caution must be taken, for example by estimating the 
memory requirement of the tasks, by grouping jobs 
based on memory, and by adapting the number of con- 
current allocated CPUs with respect to the expected 
simulation speed. In addition, in many shared high- 
performance computer facilities, walltime is limited. 
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This imposes a difficult problem because estimating 
computer time a priori is highly nontrivial, especially 
since the number of the required electronic and ionic re- 
laxations depends on how distant the initial configura- 
tion is from the unknown final equilibrium. The second 
most common cause of interruption is due to runtime 
errors of VASP. AFLOW is capable of detecting most 
of the problems and it contains many self healing fea- 
tures. This is achieved by diagnosing the error message, 
self-correcting the appropriate parameters, and restart- 
ing VASP. With AFLOW, a job can be easily restarted 
from "relaxations", "static", or "bands" steps. AFLOW's 
capability to continuously search and manage sub fold- 
ers are not limited to DFT calculation. An alien mode 
is implemented, which allows AFLOW to execute other 
tasks in a high-throughput fashion. For instance, the 
many thousands Grand Canonical Monte Carlo calcula- 
tions used in a recent surface science absorption project 



Table 1: Default value of U and J parameters given in eV ap- 
plied to f-orbitals within the GGA+U approximation included 
in AFLOW. Note that these parameters are subject to update. 



. 14711 . were directed and performed by AFLOW. In ad- 
dition, AFLOW is equipped with options to run a "pre" 
and "post" command/scripts that will be executed be- 
fore and after the main program is performed in each 
folder, respectively. This allows AFLOW to generate in- 
put files on the fly depending on the results of different 
calculations, so that ad-hoc optimization can be imple- 
mented by the users. In conclusion, the "alien" mode 
and the "pre/post" command options improves the flex- 
ibility on the recovery from a crash or an unconverged 
run in a high-throughput manner as well as increases the 
overall versatility and throughput of AFLOW. 
Implemented electronic properties. 

A typical information that one can extract from the 
band structure calculations includes the Fermi energy, 
band gap, type of the band gap, width of valence and 
conduction bands, effective mass of electron and hole, 
charge densities, band structures, total and partial den- 
sity of states, etc. A user can easily create utilities in any 
language at choice for data analysis, and use AFLOW 
in "alien" mode to execute the utilities automatically 
in each subf older s. For our purpose we have chosen 
MATLAB, which has produced all the band structure and 
orbital-projected total density of states for every shape 
the BZ as presented in the Appendix. 
LSDA+U corrections. 

It is generally known that due to a rather weak orbital- 
dependence of the DFT's exchange correlation energy, 
the strong on-site Coulomb repulsion in systems with 
narrow d- and f-bands is underaccounted. As a result, 
DFT produces bandgaps that are smaller than experi- 
mental values and some times it fails to get the cor- 
rect ground state in such systems. Based on our ex- 
perience in calculating the electronic structure of many 
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lanthanum halides, DFT incorrectly gives conduction 
band minima with 4f states instead of 5d orbitals. The 
insufficient description of strongly correlated systems 
given by DFT, can be remedied, at least partially, by 
GW 1 48] or LSDA+U corrections. Due to their large 
computational cost, GW corrections are not currently 
applicable for high-throughput searches. When needed, 
LSDA+U corrections are automatically implemented 
by AFLOW, based on the formulas developed by Du- 
radev 114911 and Liechtenstein 115011 . To the best of our 
knowledge, there are no systematic studies for the Hub- 
bard U and the screened Stoner exchange parameters 
/ across all the elements with all the possible oxida- 
tion states. Detailed analysis and determination of the 
U and / values for different compounds would be one 
of the tasks of high-throughput future research isih - In 
the mean time, we have applied the +U corrections to 
the 4f-wave functions of lanthanide compounds to get 
the correct orbitals at the conduction band minimum. 
For the systems not available in literature but related 
to our current research, Nd, Sm, and Eu, we have fit 
U and / so that the 4f levels reproduce the experimen- 
tal density of states from the x-ray photoelectron spec- 
troscopy and Bremsstrahlung isochromat spectroscopy 
(XPS-BIS) measurements D2h . Although the data is for 
metals, we are confident that the values of U and / for 
other compounds will not be very different from the fit. 
The values are listed in table [TJ 
On-line implementation: aconvasp-online. 

Users who do not need to perform high-throughput 
calculations or to create databases, can prepare standard 
unit cells input files and extract the appropriate k-points 
path by using the command version of AFLOW called 
ACONVASP or the on-line tool ACONVASP-online avail- 
able in our website (http://materials.duke.edu). 

The following protocol should be followed. Unit 
cells must first be reduced to standard primitives, then 
they should be appropriately relaxed (if needed). Be- 
fore the static run, the cells should be reduced again 
to standard primitive (symmetry and orientation might 
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have changed during the relaxation). The user should 
then perform the static run and then project the eigen- 
values along the directions which are specified in the 
"kpath" option. If the user is running AFLOW and VASP, 
the web interface can also prepare a template input file 
"aflow.in" which performs all the mentioned tasks. 



1. Appendix A 

The choice of lattice vectors implemented in AFLOW 
is given here. When the primitive lattice is the same 
as the conventional one, it is simply called "lattice". 
Variables a,b,c,a,fi,y denote the axial lengths and 
interaxial angles of the conventional lattice vectors, 
while k a ,kb,k c ,k a ,kj3,ky are those of the primitive 
reciprocal lattice vectors 01,02,03. The coordinates 
of symmetry k-points are given in fractions of bi , D2, b 3 . 

7.7. Cubic (CUB, cP) 

Lattice: 

ai = (a, 0, 0) 

a 2 = (0,a,0) 

a 3 = (0,0, a) 



Table 2: Symmetry k-points of CUB lattice. 
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Figure 1: Brillouin zone of CUB lattice. Path: T-X-M-r-R- 
X|M-R. An example of band structure using this path is given 
in Figure [26] 



7.2. Face-centered Cubic (FCC, cF) 
Conventional lattice: Primitive lattice: 
a! = ((3,0,0) a! = (0, a/2, a/2) 

a 2 = (0, a, 0) a 2 = (a/2, 0, a/2) 

a 3 = (0, 0, a) a 3 = (a/2, a/2, 0) 5 



Table 3: Symmetry k-points of FCC lattice. 
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Figure 2: Brillouin zone of FCC lattice. Path: T-X-W-K-r-L- 
U-W-L-K|U-X. An example of band structure using this path 
is given in Figure [27] 

1.3. Body-centered Cubic (BCC, cl) 
Conventional lattice: Primitive lattice: 
ai = (a, 0, 0) ai = (-a/2, a/2, a/2) 

a 2 = (0, a, 0) a 2 = (a/2, -a/2, a/2) 

a 3 = (0, 0, a) a 3 = (a/2, a/2, -a/2) 



Table 4: Symmetry k-points of BCC lattice. 
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Figure 3: Brillouin zone of BCC lattice. Path: r-H-N-r-P-H|P- 
N. An example of band structure using this path is given in 
Figure 1281 



1.4. Tetragonal (TET, tP) 

Lattice: 
a! = (a, 0, 0) 
a 2 = (0,0,0) 
a 3 = (0,0,c) 



Table 5: Symmetry k-points of TET lattice. 
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Figure 4: Brillouin zone of TET lattice. Path: r-X-MT-Z-R- 
A-Z|X-R|M-A. An example of band structure using this path 
is given in Figure [29] 



7.5. Body -centered Tetragonal (BCT, tl) 



Conventional lattice: 

a! = (a, 0, 0) 

a 2 = (0, a, 0) 

a 3 = (0,0,c) 
Variations: 
BCTi: c < a 
BCT 2 : c> a 



Primitive lattice: 
a! = (-a/2, a/2, c/2) 
a 2 = (a/2, -a/2, c/2) 
a 3 = (a/2, a/2, -c/2) 



Table 6: Symmetry k-points of BCTi lattice. 
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Figure 5: Brillouin zone of BCTi lattice. Path: T-X-M-r-Z-P- 
N-Zi-M|X-P. An example of band structure using this path is 
given in Figure [30l 



Table 7: Symmetry k-points of BCT2 lattice. 



xbi 


xb 2 


xb 3 




xbi 


xb 2 


xb 3 













r 








7 2 


X 





V2 





N 


-f 


<T 


72 


Y 


% 


% 


% 


P 


V2 


V2 


-f 


Yi 










V2 


V2 


-V2 


Z 


1 


1-7/ 


-1 


Si 










77 = 


(1 + a 2 /c 2 )/4 




f = 


a 2 /(2c 2 ) 








Figure 6: Brillouin zone of BCT 2 lattice. Path: r-X-Y-ET- 
Z-Si-N-P-Yi-Z|X-P. An example of band structure using this 
path is given in Figure [31] 
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1.6. Orthorhombic (ORC, oP) 

Ordering of the conventional lattice: a < b < c. Lattice: 
a! = (a, 0, 0) 
a 2 = (0,6,0) 
a 3 = (0, 0, c) 



Table 8: Symmetry k-points of ORC. 
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Figure 7: Brillouin zone of ORC lattice. Path: r-X-S-Y-T-Z-U- 
R-T-Z|Y-T|U-X|S-R. An example of band structure using this 
path is given in Figure [32] 



1. 7. Face-centered Orthorhombic ( ORCF, oF) 

Ordering of the conventional lattice: a < b < c. 
Conventional lattice: Primitive lattice: 
a! = (a, 0, 0) a! = (0, b/2, c/2) 

a 2 = (0, b, 0) a 2 - (a/2, 0, c/2) 

a 3 = (0,0,c) a 3 = (fl/2,&/2,0) 

Variations: 

ORCFi: I /a 2 > \/b 2 + \/c 2 
ORCF 2 : 1/a 2 <l/b 2 + 1/c 2 
ORCF 3 : 1/a 2 = 1/b 2 + 1/c 2 



Table 9: Symmetry k-points of ORCFi and ORCF 3 . 
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Figure 8: Brillouin zone of ORCFi lattice. Path: r-Y-T-Z-T-X- 
Ai-Y|T-Xi|X-A-Z|L-r. An example of band structure using 
this path is given in Figure [33] 




Figure 9: Brillouin zone of ORCF 3 lattice. Path: r-Y-T-Z-T- 
X-Ai-Y|X-A-Z|L-r. An example of band structure using this 
path is given in Figure [35] 



Table 10: Symmetry k-points of ORCF 2 . 
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Figure 10: Brillouin zone of ORCF 2 lattice. Path: T-Y-C-D-X- 
r-Z-Di-H-CICi-ZIX-HilH-YIL-r. An example of band struc- 
ture using this path is given in Figure [34] 



1.8. Body-centered Orthorhombic ( ORCI, ol) 
Ordering of the conventional lattice: a < b < c. 



Conventional lattice: 



Primitive lattice: 



ai = 


(a, 0, 0) 




ai = 


(-a/2,b/2,c/2) 




a 2 = 


(0, b, 0) 




a 2 = 


(a/2,- 


-b/2,c/2) 




a 3 = 


(0, 0, c) 




a 3 = 


(a/2,b/2, -c/2) 






Table 11: Symmetry k-points of ORCI. 
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■■ (a 2 + b 2 )/(4c 2 ) 



Table 12: Symmetry k-points of ORCC. 
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Figure 12: Brillouin zone of ORCC lattice. Path: T-X-S-R- 
A-Z-r-Y-Xi-Ai-T-Y|Z-T. An example of band structure using 
this path is given in Figure [37] 




Figure 11: Brillouin zone of ORCI lattice. Path: T-X-L-T-W- 
R-Xi-Z-r-Y-S-W|L r Y|Yi-Z. An example of band structure 
using this path is given in Figure [36] 



1.9. C-centered Orthorhombic (ORCC, oS) 

Ordering of the conventional lattice: a < b. 
Conventional lattice: Primitive lattice: 
ai = (a, 0, 0) ai = (a/2, -b/2, 0) 

a 2 = (0, b, 0) a 2 = (a/2, b/2, 0) 

a 3 = (0, 0, c) a 3 = (0, 0, c) 



1.10. Hexagonal (HEX, hP) 
Lattice: 

a! = (a/2, -(a V3)/2, 0) 
a 2 = (a/2, (a V5)/2, 0) 
a 3 = (0,0,c) 



Table 13: Symmetry k-points of HEX. 
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Figure 13: Brillouin zone of HEX lattice. Path: T-M-K-r-A-L- 
H-A|L-M|K-H. An example of band structure using this path 
is given in Figure [38] 
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7.77. Rhombohedral (RHL, hR) 
Lattice: 

ai = (a cos(of/2), -a sin(a/2), 0) 

a 2 = (acos(a/2),asm(a/2),0) 

a3 = (a cos a I cos(a/2), 0, a ^\ - cos 2 a/ cos 2 (a/2) ) 
Variations: 
RHLi: a<90° 
RHL 2 : a > 90° 




Table 14: Symmetry k-points of RHL i 
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rj = (1 + 4coso0/(2 + 4 cos a) 
v= 3 / 4 -77/2 




Figure 14: Brillouin zone of RHLi lattice. Path: T-L-Bi|B-Z-r- 
X|Q-F-Pi-Z|L-R An example of band structure using this path 
is given in Figure [39] 



Table 15: Symmetry k-points of RHL2. 
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Figure 15: Brillouin zone of RHL 2 lattice. Path: T-P-Z-Q-r- 
F-P1-Q1-L-Z. An example of band structure using this path is 
given in Figure |40l 



7.72. Monoclinic (MCL, mP) 

Ordering of the lattice: a, b < c, a < 90°, J3 = y = 90° 

Lattice: 

ai = (a, 0, 0) 

a 2 = (0, b, 0) 

a 3 = (0, c cos a, c sin a) 

Table 16: Symmetry k-points of MCL. 
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77 = (1 - b cos a/c) I (2 sin a) 
v = V2 - cos a lb 




Figure 16: Brillouin zone of MCL lattice. Path: T-Y-H-C-E- 
Mi-A-X-Hi|M-D-Z|Y-D. An example of band structure using 
this path is given in Figure I4T1 
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1.13. C-centered Monoclinic (MCLC, mS) 



Ordering of the conventional lattice: a,b < c,a < 

90°,/? = r = 90°. 
Conventional lattice: Primitive lattice: 
ai = (a, 0,0) ai = (a/2,b/2,0) 

a 2 = (0, b, 0) a 2 = (-a/2, b/2, 0) 

a3 = (0, c cos a, c sin a) a3 = (0, c cos a, c sin a) 

Variations: 
MCLCi: k 7 > 90° 
MCLC 2 : iky = 90° 

MCLC 3 : k y < 90°,bcosa/c + b 2 sin 2 a/ a 2 < 1 
MCLC 4 : k y < 90°,bcosa/c + b 2 sin 2 a/a 2 = 1 
MCLC 5 : ky < 90°,bcosa/c + b 2 sin 2 a/a 2 > 1 




Figure 18: Brillouin zone of MCLC 2 lattice. Note that Y is 
equivalent to X. Path: r-Y-F-L-I|I r Z-Fi|N-r-M. An example 
of band structure using this path is given in Figure [43] 



Table 17: Symmetry k-points of MCLCi and MCLC 2 . 
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f = (2-Z?cosa/c)/(4sin 2 a) 
7] = \ + 2gccosa/b 
ij/ = % - a 2 /(4b 2 sin 2 a) 
cp = if/ + ( 3 / 4 - cos ar/c 



Table 18: Symmetry k-points of MCLC 3 and MCLC 4 . 
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ja = (1 + b 2 /a 2 )/4 
S = be cos a /(2a 2 ) 

£ = /j, - V 4 + (1 - /?cosa/c)/(4sin 2 a) 
7/ = V 2 + 2£ccosa/b 
= 1 + f - 2/i 
if/ = 7] -26 




I-Fi|Hi-Yi-X-r-N|M-r. An example of band structure using 
Figure 17: Brillouin zone of MCLCi lattice. Path: T-Y-F-L- this path is given in Figure [44] 
I|Ii-Z-Fi|Y-Xi|X-r-N|M-r. An example of band structure us- 
ing this path is given in Figure [42] 




Figure 20: Brillouin zone of MCLC 4 lattice. Note that I is 
equivalent to F. Path: T-Y-F-H-Z-IIHrYrX-r-NIM-r. An ex- 
ample of band structure using this path is given in Figure [45] 



Table 19: Symmetry k-points of MCLC 5 
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7.74. Triclinic (TRI, aP) 

Lattice: 

ai = (a, 0, 0) 

a 2 = (b cos 7, b sin 7, 0) 

a 3 = (c cos /?, ^ [cos of - cos J3 cos 7] , 

^- A / sin 2 7 - cos 2 a - cos 2 /? + 2 cos a cos/3 cos 7 ) 

siny y 

Variations: 



TRIla 

TRIi, 
TRI 2 « 
TRI 2 , 



k a >90 o ,k fi >90°,ky >90° 
^ < 90°, < 90°, £ r < 90° 
Jfca > 90°, > 90°, £ r = 90° 
k a < 90°, kj3< 90°, ky = 90° 



Table 20: Symmetry k-points of TRI ia and TRI 2fl . 
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Figure 21: Brillouin zone of MCLC 5 lattice. Path: T-Y-F-L- 
I|Ii-Z-H-Fi|Hi-Yi-X-r-N|M-r. An example of band structure 
using this path is given in Figure [46] 



Figure 22: Brillouin zone of TRIi a lattice. Path: X-r-Y|L-T- 
Z|N-r-M|R-r. An example of band structure using this path is 
given in Figure [47] 
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Figure 23: Brillouin zone of TRI 2a lattice. Path: X-r~Y|L-r- 

z|N-r-M|R-r. 



2. Appendix B 



Table 21: Symmetry k-points of TRI^ and TRI 2 ^. 
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Figure 24: Brillouin zone of TRI^ lattice. Path: X-r-Y|L-r- 
Z|N-r-M|R~r. An example of band structure using this path is 
given in Figure [48] 



Selected examples of band structure for each shape 
of Brillouin zone are presented in this section. The 
Fermi energy is shifted to the valence band maximum at 
zero. In each figure, the orbital-projected total density 
of states N(E) are plotted in the right panel in logarith- 
mic scale. 



PaSr^Sr, ICSD 27047 (CUB) 




" 10 r X M r R X|M R log(N(E)) 

Figure 26: Band structure of Sr(Sn0 3 ) in CUB lattice. 



Cd 1 S 1 ICSD 29278 (FCC) 




" 10 r xwk r l uw lkiux iog(N(E» 

Figure 27: Band structure of CdS in FCC lattice. 
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2 Sn^ ICSD9163 (TET) 




l0 r x m r z r A Z|X R|M a iog(N(E)) 

Figure 29: Band structure of Sn0 2 in TET lattice. 



O^Zn,, ICSD9147 (BC^) 




'°r x M r z P N ZiM|x P iog(N(E)) 

Figure 30: Band structure of Zn 2 (Si0 4 ) in BCTi lattice. 



Ca^n^ ICSD 52390 (BCT 2 ) 




10 r X Y E r Z £, N P Yi Z|X P log(N(E)) 

Figure 31: Band structure of Caln 2 4 in BCT 2 lattice. 



Bi t F 3 ICSD 1269 (ORC) 




10 r X S Y r Z U R T Z|Y T|U X|S R log(N(E)) 

Figure 32: Band structure of BiF 3 in ORC lattice. 



CdjC^Si, ICSD 23943 (ORCF,) 




" 1 °r Y T z r X A YITX^X A Z|L r log(N(E)) 

Figure 33: Band structure of Cd 2 (Si0 4 ) in ORCFi lattice. 



F 1 ICSD 30268 (ORCF 2 ) 




" 1 °r YCDX r Z D] H C]Ci Z|XH-||H Y|L r log(N(E)) 

Figure 34: Band structure of TIF in ORCF 2 lattice. 



F 1 ln 1 1 ICSD 2521 (ORCF 3 ) 




-8 - ; ; ; ; ; ; ; j < 

10 r Y T z r X Ai Y|X A Z|L r log(N(E)) 

Figure 35: Band structure of InOF in ORCF 3 lattice. 



Se^ ICSD 24592 (ORCI) 




-8- ; ; ; ; ill ! ! ! ! ! 

'°r XLTW RX!Z r YS W|L, Y|Y z iog(N(E)) 

Figure 36: Band structure of SiSe 2 in ORCI lattice. 
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1^ ICSD 26761 (ORCC) 




Figure 41: Band structure of Bi(B0 3 ) in MCL lattice. 



Bi^ ICSD 20676 (HEX) 




Figure 38: Band structure of Bil 3 in HEX lattice. r Y FL 11,1 z F i |YX ^ |X r N|M r l09(N(E)) 

Figure 42: Band structure of BaNa(B0 3 ) in MCLCi lattice. 



AI 1 Bi 1 3 ICSD 171708 (RHL^) 




Figure 39: Band structure of BiA10 3 in RHLi lattice. 



Ba 1 F 6 Ge 1 ICSD 26614 (RHL 2 ) 




'°r PZ Q r F ^ Q 1 L Z log(N(E)) 

Figure 40: Band structure of Ba(GeF 6 ) in RHL 2 lattice. 



La 2 6 Si ICSD 66823 (MCLC 2 ) 




" 10 r Y F L III, Z RJN r M log(N(E)) 

Figure 43: Band structure of (LaO) 2 (S0 4 ) in MCLC 2 lattice. 



Br 3 ln 1 ICSD 65198 (MCLC 3 ) 




-6 
-8 



10 r Y F H Z I F 1 |H 1 Y, X r N|M r log(N(E)) 

Figure 44: Band structure of InBr 3 in MCLC 3 lattice. 
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Na 4 9 Si 3 Sr 1 ICSD 33943 (MCLC 4 ) 




" 10 r Y F HZ l|H-, V, X r N|M r log(N(E)) 

Figure 45: Band structure of Na 4 Sr(Si0 3 ) 3 in MCLC 4 lattice. 

Cd 1 Na 2 S 4 Sn 1 ICSD 281233 (MCLC 5 ) 




-8- !!!!!! ! ! ! ! ! 

" 10 r Y F L l|h Z H F 1 |H 1 Y, x r N|M r log(N(E)) 

Figure 46: Band structure of Na 2 (CdSnS 4 ) in MCLC 5 lattice. 



B 1 Cd 1 Li 1 3 ICSD 200615 (TRI 1a ) 




-8 



10 x r Y|L r Z|N r M|R r iog(N(E)) 

Figure 47: Band structure of LiCd(B0 3 ) in TRI ia lattice. 



B 2 Cd 2 5 ICSD 281357 (TRI 1b ) 




'°x r Y|L r Z|N r M|R r iog(N(E)) 

Figure 48: Band structure of Cd 2 (B 2 5 ) in TRI^ lattice. 
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